CASPON platform technology: Ultrafast circularly permuted caspase-2 cleaves tagged fusion proteins before all 20 natural amino acids at the N-terminus.

Fusion protein technologies improve the expression and purification of recombinant proteins, but the removal of the tags involved requires specific proteases. The circularly permuted caspase-2 (cpCasp2) with its specific cleavage site, efficiently generates the untagged protein. While cleavage with cpCasp2 is possible before all 20 proteinogenic amino acids, cleavage before valine, leucine, isoleucine, aspartate and glutamate suffers from slow, and before proline extremely slow, turnover. To make the platform fusion protein process even more general such that any protein with an authentic N-terminus can be produced with high efficiency, the bacterial selection system PROFICS (PRotease Optimization via Fusion-Inhibited Carbamoyltransferase-based Selection) was used to evolve cpCasp2 into a variant with a catalytic turnover two orders of magnitude higher and the ability to cleave before any amino acid. The high specificity and the stability of the original circularly permuted protease was fully retained in this mutant, while the high manufacturability was mostly retained, albeit with decreased soluble titer. Four point-mutations are responsible for this change in activity, two of which are located in or near the binding pocket of the active site. This variant was named CASPON enzyme and is a major component of the CASPase-based fusiON (CASPON) platform technology. Applicability for the production of recombinant proteins was demonstrated by enzymatic removal of the CASPON tag from five model proteins. The CASPON tag enables high soluble expressions, affinity purification and good accessibility for cleavage. The five industry-relevant proteins of interest were FGF2, TNF, GH, GCSF and PTH.

BuRNN: Buffer Region Neural Network Approach for Polarizable-Embedding Neural Network/Molecular Mechanics Simulations.

Hybrid quantum mechanics/molecular mechanics (QM/MM) simulations have advanced the field of computational chemistry tremendously. However, they require the partitioning of a system into two different regions that are treated at different levels of theory, which can cause artifacts at the interface. Furthermore, they are still limited by high computational costs of quantum chemical calculations. In this work, we develop the buffer region neural network (BuRNN), an alternative approach to existing QM/MM schemes, which introduces a buffer region that experiences full electronic polarization by the inner QM region to minimize artifacts. The interactions between the QM and the buffer region are described by deep neural networks (NNs), which leads to the high computational efficiency of this hybrid NN/MM scheme while retaining quantum chemical accuracy. We demonstrate the BuRNN approach by performing NN/MM simulations of the hexa-aqua iron complex.

Reaction intermediate rotation during the decarboxylation of coproheme to heme b in C. diphtheriae.

Monoderm bacteria utilize coproheme decarboxylases (ChdCs) to generate heme b by a stepwise decarboxylation of two propionate groups of iron coproporphyrin III (coproheme), forming two vinyl groups. This work focuses on actinobacterial ChdC from Corynebacterium diphtheriae (CdChdC) to elucidate the hydrogen peroxide-mediated decarboxylation of coproheme via monovinyl monopropionyl deuteroheme (MMD) to heme b, with the principal aim being to understand the reorientation mechanism of MMD during turnover. Wild-type CdChdC and variants, namely H118A, H118F, and A207E, were studied by resonance Raman and ultraviolet-visible spectroscopy, mass spectrometry, and molecular dynamics simulations. As actinobacterial ChdCs use a histidine (H118) as a distal base, we studied the H118A and H118F variants to elucidate the effect of 1) the elimination of the proton acceptor and 2) steric constraints within the active site. The A207E variant mimics the proximal H-bonding network found in chlorite dismutases. This mutation potentially increases the rigidity of the proximal site and might impair the rotation of the reaction intermediate MMD. We found that both wild-type CdChdC and the variant H118A convert coproheme mainly to heme b upon titration with H2O2. Interestingly, the variant A207E mostly accumulates MMD along with small amounts of heme b, whereas H118F is unable to produce heme b and accumulates only MMD. Together with molecular dynamics simulations, the spectroscopic data provide insight into the reaction mechanism and the mode of reorientation of MMD, i.e., a rotation in the active site versus a release and rebinding.

PROFICS: A bacterial selection system for directed evolution of proteases.

Proteases serve as important tools in biotechnology and as valuable drugs or drug targets. Efficient protein engineering methods to study and modulate protease properties are thus of great interest for a plethora of applications. We established PROFICS (PRotease Optimization via Fusion-Inhibited Carbamoyltransferase-based Selection), a bacterial selection system, which enables the optimization of proteases for biotechnology, therapeutics or diagnosis in a simple overnight process. During the PROFICS process, proteases are selected for their ability to specifically cut a tag from a reporter enzyme and leave a native N-terminus. Precise and efficient cleavage after the recognition sequence reverses the phenotype of an Escherichia coli knockout strain deficient in an essential enzyme of pyrimidine synthesis. A toolbox was generated to select for proteases with different preferences for P1' residues (the residue immediately following the cleavage site). The functionality of PROFICS is demonstrated with viral proteases and human caspase-2. PROFICS improved caspase-2 activity up to 25-fold after only one round of mutation and selection. Additionally, we found a significantly improved tolerance for all P1' residues caused by a mutation in a substrate interaction site. We showed that this improved activity enables cells containing the new variant to outgrow cells containing all other mutants, facilitating its straightforward selection. Apart from optimizing enzymatic activity and P1' tolerance, PROFICS can be used to reprogram specificities, erase off-target activity, optimize expression via tags/codon usage, or even to screen for potential drug-resistance-conferring mutations in therapeutic targets such as viral proteases in an unbiased manner.

Charge-Changing Perturbations and Path Sampling via Classical Molecular Dynamic Simulations of Simple Guest-Host Systems.

Currently, two different methods dominate the field of biomolecular free-energy calculations for the prediction of binding affinities. Pathway methods are frequently used for large ligands that bind on the surface of a host, such as protein-protein complexes. Alchemical methods, on the other hand, are preferably applied for small ligands that bind to deeply buried binding sites. The latter methods are also widely known to be heavily artifacted by the representation of electrostatic energies in periodic simulation boxes, in particular, when net-charge changes are involved. Different methods have been described to deal with these artifacts, including postsimulation correction schemes and instantaneous correction schemes (e.g., co-alchemical perturbation of ions). Here, we use very simple test systems to show that instantaneous correction schemes with no change in the system net charge lower the artifacts but do not eliminate them. Furthermore, we show that free energies from pathway methods suffer from the same artifacts.

Actinobacterial Coproheme Decarboxylases Use Histidine as a Distal Base to Promote Compound I Formation.

Coproheme decarboxylases (ChdCs) catalyze the final step in heme b biosynthesis of monoderm and some diderm bacteria. In this reaction, coproheme is converted to heme b via monovinyl monopropionate deuteroheme (MMD) in two consecutive decarboxylation steps. In Firmicutes decarboxylation of propionates 2 and 4 of coproheme depend on hydrogen peroxide and the presence of a catalytic tyrosine. Here we demonstrate that ChdCs from Actinobacteria are unique in using a histidine (H118 in ChdC from Corynebacterium diphtheriae, CdChdC) as a distal base in addition to the redox-active tyrosine (Y135). We present the X-ray crystal structures of coproheme-CdChdC and MMD-CdChdC, which clearly show (i) differences in the active site architecture between Firmicutes and Actinobacteria and (ii) rotation of the redox-active reaction intermediate (MMD) after formation of the vinyl group at position 2. Distal H118 is shown to catalyze the heterolytic cleavage of hydrogen peroxide (k app = (4.90 ± 1.25) × 104 M-1 s-1). The resulting Compound I is rapidly converted to a catalytically active Compound I* (oxoiron(IV) Y135•) that initiates the radical decarboxylation reactions. As a consequence of the more efficient Compound I formation, actinobacterial ChdCs exhibit a higher catalytic efficiency in comparison to representatives from Firmicutes. On the basis of the kinetic data of wild-type CdChdC and the variants H118A, Y135A, and H118A/Y135A together with high-resolution crystal structures and molecular dynamics simulations, we present a molecular mechanism for the hydrogen peroxide dependent conversion of coproheme via MMD to heme b and discuss differences between ChdCs from Actinobacteria and Firmicutes.

Correcting electrostatic artifacts due to net-charge changes in the calculation of ligand binding free energies.

Alchemically derived free energies are artifacted when the perturbed moiety has a nonzero net charge. The source of the artifacts lies in the effective treatment of the electrostatic interactions within and between the perturbed atoms and remaining (partial) charges in the simulated system. To treat the electrostatic interactions effectively, lattice-summation (LS) methods or cutoff schemes in combination with a reaction-field contribution are usually employed. Both methods render the charging component of the calculated free energies sensitive to essential parameters of the system like the cutoff radius or the box side lengths. Here, we discuss the results of three previously published studies of ligand binding. These studies presented estimates of binding free energies that were artifacted due to the charged nature of the ligands. We show that the size of the artifacts can be efficiently calculated and raw simulation data can be corrected. We compare the corrected results with experimental estimates and nonartifacted estimates from path-sampling methods. Although the employed correction scheme involves computationally demanding continuum-electrostatics calculations, we show that the correction estimate can be deduced from a small sample of configurations rather than from the entire ensemble. This observation makes the calculations of correction terms feasible for complex biological systems. To show the general applicability of the proposed procedure, we also present results where the correction scheme was used to correct independent free energies obtained from simulations employing a cutoff scheme or LS electrostatics. In this work, we give practical guidelines on how to apply the appropriate corrections easily.